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The journey through the nontrivial band 
topology beyond the conventional band structure 
has resulted in the recent discovery of three- 
dimensional Dirac semimetal phase in Na 3 Bi 
and Cd 3 As 2 [1-7]. The bulk state of which is 
semi-metallic obeying linear energy dispersion, 
while the surface state is topology protected 
Fermi arc [1, 2, 8, 9], Due to the unique band 
topology, they show different exotic electronic 
properties of both fundamental and technological 
interest. From electrical and thermal trans¬ 
port measurements, we have demonstrated a 
remarkable violation of Wiedemann-Franz law 
(WFL) under application of magmatic field in 
Cd 3 As 2 and the violation becomes more and 
more drastic with increasing magnetic field 
strength. Whereas the validity of WFL is the key 
feature of Landau Fermi-liquid theory in metal, 
the notion of quasiparticles is the building block 
to this theory. This implies that the fundamental 
concept of Landau quasiparticle no longer holds 
in Cd 3 As 2 in presence of magnetic field. The 
continuous break down of Landau quasiparticle 
framework with field introduces a concept of 
field induced quantum critical point (QCP) as in 
the case of heavy fermion compounds YbRh 2 Si 2 
[10], Sr 3 Ru 2 07 [11], etc. 

Dirac fermionic excitation and topology introduce 
different exotic electronic properties of solids which 
are now the subject of considerable research interest 
in condensed matter physics, both from basic research 
and application point of view. This unique excitation 
was first discovered in two-dimensional graphene [12] 
followed by topological insulator [13, 14] and the parent 
compounds of iron-based superconductors [15]. Un¬ 
like conventional fermionic excitations in metals, the 
quasiparticle excitation in these systems follows a linear 
energy dispersion. Recently, the three-dimensional 
(3D) Dirac semimetal state has been theoretically 
predicted, in which the semimetallic bulk is the 3D 
analog of graphene [1, 2, 8, 9]. The 3D Dirac points, 
where the two Weyl points overlap in momentum 
space, are protected by the crystal symmetry [1, 2, 
8, 9]. Breaking of time reversal symmetry by exter¬ 
nal magnetic field ( B) rearranges the Fermi surface [1, 2]. 

Following theoretical prediction [1, 2], 3D Dirac 
semimetal phase has been discovered experimentally in 
Na 3 Bi and Cd 3 As 2 [3-7]. So far, the major attentions 


have been focused on the electronic properties of Cd 3 As 2 
for understanding the nature of band structure [16-19]. 
It has been observed that the Fermi surface of Cd 3 As 2 
undergoes a change with the application of magnetic field 
[18]. In case of heavy fermion compounds like YbRh 2 Si 2 , 
the Fermi surface is strongly affected by B [10, 20]. 
Above a critical field ( B c ), the paramagnetic phase is a 
heavy Fermi liquid with large Fermi surface due to the 
Kondo effect and the antiferromagnetic phase below B c 
is also a mass enhanced Fermi liquid with small Fermi 
surface [20]. At B=B C , however, a non-Fermi-liquid 
ground state emerges in the field-induced QCP as T 
approaches towards zero [10]. This motivated us to 
investigate whether the concept of Landau Fermi-liquid 
holds across the Fermi surface reconstruction in presence 
of magnetic field in Cd 3 As 2 , where instead of strong 
electron-electron correlation, the spin-orbit coupling 
plays a crucial role [1]. 


Here, we focus on the electrical and thermal transport 
properties of Cd 3 As 2 single crystal in magnetic field. 
Figure 1 displays the normalized magnetoresistivity (p) 
at different temperatures, p increases almost linearly 
with increasing held. Even at 300 K and 9 T magnetic 
held, p shows no sign of saturation, and the resistivity 
ratio p{B)/p{ 0) is~3. With decreasing temperature, 
the value of p(B)/p( 0) increases rapidly and at 10 K, 
it is ~17.5 at 9 T. Inset of Fig. 1 shows temperature 
dependence of p. p does not follow T 2 and linear T 
dependence at low temperature and high temperature 
respectively as in the case of normal metal. 


In Fig. 2 (a), n has been plotted as a function of B 
at some selected temperatures between 20 and 300 K. 
k decreases very rapidly with increasing B and tends to 
saturate at high B which is expected to be the lattice 
part of thermal conductivity (kl)- This saturating be¬ 
haviour is more prominent below 50 K. The inset of Fig. 
2 (a) shows temperature dependence of k. The sharp in¬ 
crease in k below 50 K also suggests the smaller propor¬ 
tion of electronic thermal conductivity (re e ) to the total 
at low temperature. We have fitted the field dependence 
of thermal conductivity with 


k = k l 


K e 

1 + p\B 2 


(1) 


to separate n e from total thermal conductivity k, where 
the B—^oo limit will give the value of kl and px is 
the thermal mobility [21-24]. At low temperature, the 
extracted value of kl is close to the value of k at 9 T. 




FIG. 1: (color online) Resistivity (p) of CCI 3 AS 2 single crystal 
in presence of magnetic field at some selected temperatures. 
Inset shows the temperature dependence of p in absence of 
field. 


As an example, at 20 K, we have deduced kl— 3.549 
W K -1 m -1 which is close to the observed value of re 
at 9 T (3.55 W K -1 m -1 ). If one considers the 1st 
order expansion of the above equation for 1 (for 

very small B), one finds K-KL=K e (B)=K e (l-fiT 2 B 2 ). 
For a 3D Dirac semimetal, the quadratic B of re e has 
been theoretically predicted by considering an energy 
independent scattering time (r), which is given by 
K e {B) oc [1 — (w c t) 2 ], where w c ( =eB/Ti ) is the cyclotron 
frequency of the electron orbit [25]. This implies that the 
above assumption of energy independence of scattering 
time is not valid at high field. 

Below 20 K, instead of showing saturation, re decreases 
almost linearly up to 9 T and this effect enhances in 
the phonon dominated region of re, as shown in Fig. 
2 (b). Now, equation (1) can not be used to separate 
re l from re. We have measured magnetostriction at low 
temperature to check whether B induces any structural 
change. But we have not observed any noticeable 
changes in lattice parameters. The specific heat does not 
show any B dependence which implies that the phonon 
density of state is not affected by B. So, we propose 
that the observed phenomenon is due to the reduction in 
lattice thermal conductivity in magnetic field. Recently, 
the lattice thermal conductivity of diamagnetic semicon¬ 
ductor InSb is reported to exhibit B dependence [26]. 
In that work, the authors have argued that B creates 
a spatial gradient in diamagnetic moment around the 
displaced atoms in phonon vibration which exerts an 
anharmonic magnetic force on the displaced atom and 
as a result, a reduction in lattice thermal conductivity. 
We have estimated about 10% reduction in kl at 6 K 
and 9 T which is comparable to that for InSb [26]. 


t - 1 -1-'-1-'-r 
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FIG. 2: (color online) Magnetic field dependence of thermal 
conductivity (a) From 20 to 300K at some representative tem¬ 
peratures up to 9 T. Solid lines are the fitted to the exper¬ 
imental data. Inset shows the temperature dependence of 
thermal conductivity at zero held, (b) Below 20 K up to 9 
T. It has been fitted by equation (1) at 15 and 10 K over 
a smaller region (solid lines) relative to higher temperatures 
data. 


To calculate re e /crT, i.e, the Lorentz number (L), we 
have extracted re e from re using equation 1. Figure 3 (a) 
displays the temperature variation of L in absence of 
magnetic field. L over most of the temperature range 
is only about 70% of the ideal value Lo (2.44xl0 -8 
WDK -2 ). L shows a shallow dip around 100 K and 
approaches rapidly to Lo at low temperature. Below 
10 K, it is difficult to calculate the precise value of 
L because kl is affected by B. The temperature 
dependence of L is qualitatively similar to earlier report 
[22], where the B dependence of re was measured only up 
to 2 T to separate re e using equation 1 and to calculate L. 

Figure 3 (b) shows n e /aT ratio under application 
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FIG. 3: (color online) (a) Temperature dependence of Lorentz 
number (L), K e /aT, at zero field, (b) Reduction of Ke/aT 
ratio under application of magnetic field upto 9 T at some 
constant temperatures. 


of B up to 9 T at some selected temperatures in the 
range 10 to 300 K. K e /aT decreases very rapidly with 
increasing B. At 300 K, it decreases from 1.95 to 0.37 
as B increases from zero to 9 T. At low temperatures, 
however, L reduces more rapidly with increasing field 
strength. For example, at 20 K, n e /aT drops from its 
zero Held value 2.39 to ^0.006, only at 3 T. There may 
be some error in the calculated value of L at 10 K due 
to smaller proportion of n e and the effect of B on kl 
but the qualitative behavior of L(B ) remains same. 
At high fields and low temperature, L is found to be 
few orders of magnitude smaller than Lo- This clearly 
demonstrates a drastic violation of the Wiedemann- 
Franz law under application of magnetic field in Dirac 
semimetal CdaAs 2 . The violation of this universal law 
can be observed from the non-identical B dependence of 
k and a. p(B) does not follow the predicted quadratic 
B dependence, p(B)/p( 0) = 1 + PpB 2 , but linear in B 



FIG. 4: (color online) Figure shows the effect of magnetic field 
on Seebeck coefficient ( S ) upto 9 T at some representative 
temperatures between 20 to 300 K. Right side of downward 
arrow represents the value of m close to 1. 


[23, 24]. As a result of such unusual linear B depen¬ 
dence of p 1 the ratio K e (B)/a(B) fails to remain constant. 

In order to understand the effect of B on the Fermi 
surface, we have measured the Seebeck coefficient ( S ). 
S increases monotonically with B and tends to saturate 
at high B as shown in Fig. 4. Also with decreasing 
temperature, the value of S decreases monotonically and 
tends to saturate at a relatively lower field strength. The 
well known Mott’s semiclassical formula of thermoelectric 
power, S = dC g^ | e =/x; h as been used in graphene 

to probe the relaxation process [27]. Here d(e) is energy 
dependent conductivity, p, is the chemical potential and 
p=Ep for T«Tp. Using the similar algorithm in case 
of ”3D analog of graphene” and scattering time (t)oc £ m , 
the above equation simplifies to, S = n ^ e B + 2). 

From the linear T dependence of S , we have deduced 
scattering exponent ro~0.15 for the relaxation process in 
CdsAs 2 single crystal with Fermi energy ~270 meV [28]. 
This linear T dependence of S is robust under application 
of field (see Supplementary Information). Now, consider¬ 
ing the above equation and Fig. 4, it can be shown that 
to increases monotonically with B and tends to saturate 
at 1 at high field. Also, it is apparent from Fig. 4 that 
m reaches close to 1 at lower B with decreasing tempera¬ 
ture. The value of to determines the nature of relaxation 
process of the charge carrier. For examples, m =2 for 
screened charged impurity scattering, m =-2 for short- 
range disorder [25] and for unscreened charged impurity 
scattering to is 1 [27], etc. The strong B dependence of 
to can be ascribed to field induced change of Fermi sur¬ 
face, similar to the earlier report [18]. 

With increasing magnetic field strength, the scattering 
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FIG. 5: (color online) A schematic representation. Here, A 
representing the region where the value of Lorentz number (L) 
is larger than in region B. Deeper shaded region represents 
the value of L closer to Lo and at (0,0), L=Lo- Also it should 
be noted that there is no sharp boundary between region A 
and B. System continuously evolves from region A to region 
B with the variation of magnetic field and temperature. 


exponents m tends to saturate at 1 while the Lorentz 
number L falls drastically from its ideal value Lq. This 
implies a continuous violation of Wiedemann-Franz 
law (WFL), i.e, the Landau quasiparticle no longer 
holds under application of magnetic field. However, 
the present and earlier studies [22] clearly show that L 
approaches Lo at low temperature and zero field. Based 
on our experimental observation, we have constructed a 
schematic diagram in Fig. 5. A represents the shaded 
region, where L is closer to Lq than the value in region B 
and at the origin (0,0), L=Lq. The system continuously 
evolves from A to B and vice versa with the variation 
of m(B). Now, by tuning m from ~ 0 to 1 through 
magnetic field, one observes a continuous breakdown of 
Landau Fermi liquid behaviour. Thus, m=l looks like a 
quantum critical point (QCP) at T tends to zero limit. 
As if, we have been able to achieve the one side of the 
point. The other side can not be accessed by tuning B 
as m saturates at 1. Similar kind of behaviour has been 
observed in heavy fermion compound like YbRli 2 Si 2 
[10], in which B induces a QCP between paramagnetic 
heavy Fermi-liquid state and antiferromagnetic mass 
enhanced Fermi-liquid state, where non-Fermi-Liquid 
ground state emerges. The present study introduces 
the concept of field-induced QCP in Dirac semimetal 
CdsAs 2 , where spin-orbit coupling is more important 
than electron-electron interaction. A qualitatively 
similar kind of behaviour has been predicted using the 
renormalization group analysis [29]. In that work, the 
authors have claimed that a quantum phase transition 
in 3D Dirac semimetals occurs due to the fascinating 
interplay of Coulomb interaction and disorder. The 


non-Fermi liquid ground state can be achieved by in¬ 
creasing the strength of disorder. This theoretical study 
is consistence with our observation. It seems that B is 
playing the role of disorder in the present study. As an 
another possible explanation, the field-induced violation 
of conventional WFL may be an interesting property 
of topological Fermi-liquid [30] (Landau Fermi-liquid 
including Berry’s phase in presence of broken inversion 
or time-reversal symmetry), which is unexplored till 
now. A quadratic B dependence of Lorentz number has 
been theoretically predicted for Weyl semimetal [31] 
(which is also a topological Fermi-liquid system), due 
to the chiral anomaly. Further theoretical studies are 
necessary to understand the possible origin of such a 
strong violation of WFL. 

Methods 

Single crystals of Cd 3 As 2 were synthesized by chemical 
vapor transport technique [28]. Powder X-ray diffraction 
of crashed single crystals shows that these crystals have 
hi/acd space group and contain no impurity phases [28]. 
The CdaAs 2 single crystal was cut and polished to a bar 
shape, with dimension^3x2x0.6x mm 3 for transport 
measurements. The resistivity, thermal conductivity 
and thermoelectric power (Seebeck coefficient) measure¬ 
ments on CdsAs 2 single crystals were done by four-probe 
technique in a physical property measurement system 
(Quantum Design) upto 9 T. Though several single 
crystals have been studied, we present the data for a 
single crystal as a representative. Qualitative similar 
behavior has been observed for other crystals. Creating 
a temperature gradient along the current direction, the 
measurements were performed by applying magnetic 
field perpendicular to them. 
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Supplementary information for ’’Mag¬ 
netic field induced drastic violation of 
Wiedemann-Franz law in Dirac semimetal 
Cd 3 As 2 ” 

I. TEMPERATURE DEPENDENCE OF 
SEEBECK COEFFICIENT (S) UNDER 
APPLICATION MAGNETIC FIELD (B). 


Figure 6 shows linear T dependence of S 
both in presence and absence of external mag¬ 
netic held ( B ) upto 350K. This implies Mott’s 
semiclassical formula of thermoelectric power 
also holds in presence of B. So the expression, 
S = + 2), can be used to calculate 

the B dependence of scattering exponent (m). 



FIG. 6: (Color online) Temperature dependence of S at OT, 
5T and 9T magnetic field. 
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